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SUMMARY 


A theory for obtaining approximate solutions to nonlinear problems whose 
exact solutions require the use of large computational procedures is described. 
The technique represents in some respects a generalization of the method of 
base and comparison solutions for flows depending on a parameter. For the 
generalized problem, the input variable is no longer a parameter but a function 
that is incremented over its entire domain. After performing calculations for 
a base configuration and a small number of variations of it, solutions for a 
large class of configurations can be obtained by forming linear combinations of 
the solution increments. For a restricted class of problems, approximate solu- 
tions can be obtained for general variations of a base configuration by using 
a function-space derivative estimate obtained frati a base solution and a single 
variation. 


INTRODUCTICW 

Many important aerospace-related problems require the solution of one or 
more nonlinear differential equations. Although analytic solutions are rare, 
many such problems have been solved by means of numerical computerized proce- 
dures. Often such procedures are long and expensive, when large computing 
times and storage are required, and they are cumbersome to operate when many 
input variables and parameters are required. 

An alternate to using such programs is to seek approximate solutions by 
linearizing the differential equations. Such a linearization procedure consists 
of considering the problem as a perturbation of a "base" problem for which the 
solution is known. However, even the linear equations often defy solution in 
analytic form unless the base problem is almost trivially simple, as for example, 
an undisturbed free-stream flow. 

Clearly the usage of large complex computational algorithms could be 
reduced if there were a way to combine a program for computing a base solution, 
with a fast, simple procedure for computing the effects of perturbation about 
that solution. Such a method has been described in reference 1 for flows that 
depend on some parameter such as Mach number, angle of attack, or thickness 
ratio, when one of these parameters is perturbed from its base value. 

The present analysis represents, in seme sense, an extension of the method 
of reference 1 to problems that involve the perturbation, not of a parameter, 
but of an entire function. Thus, the method applies to problems for which the 
input is a function, such as one prescribing a boundary shape; and the output 
is also a function, such as the pressure distribution over the boundary. The 
method could also be used to determine variations in the ground noise pattern 
of an airplane due to variations in its trajectory. 





Advantages of the method are that it is applicable to a large class of 
complex problems, and that it does not require tedious analytic procedures such 
as solving inner and outer problems, matching, etc. Its primary limitation is 
that it requires a means of generating base solutions. Furthermore, it is, in 
general, not directly applicable to problems containing jump discontinuities, 
although it can be modified to handle such problems by methods similar to those 
of reference 2. 


SYMBOLS 

a constant coefficient 

Cp pressure coefficient 

c local sound speed 

F solution (output) function 

f independent (input) function 

G Gateaux derivative 

h continuous function of small maximum ordinate representing an 

increment to f 

L denotes a linear operator 

M Mach number 

p parameter 

Q denotes a nonlinear operator 

r nondimens ional body radius (in terms of body length) 

x,y nondimensional distance in free-stream direction and perpendicular 

to free-stream direction, respectively (in terms of body length) 

X shifted coordinate 

a angle of attack 

A before a variable denotes an increment in the variable 

6 before a function denotes a perturbation in the function 

e small parameter 

<j) velocity potential 
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Subscripts: 

i represents ith quantity in a sequence 

o denotes base value 

1,2,3 denotes quantities associated respectively with the first, second, 
and third variations 

x,y denotes derivatives with respect to x and y 

Superscript: 

' Gateaux derivative of an operator 


DEVELOPMENT OF THEORY 

Background Theory on Nonlinear Problems Depending on a Parameter 

Reference 1 treats the problem of reducing computer time by combining 
calculations which have a full nonlinear numerical solution with calculations 
obtained by a simple linear interpolation scheme. If the solution depends on 
a parameter p (e.g., Mach number, angle of attack, thickness ratio), then the 
effect on the solution F of an incremental change in p is 

9f 

Af = — Ap (1 ) 

9p 


provided that F does not change discontinuously with p. The latter situation 
would occur, for example, in certain transonic flows involving shock waves. The 
problem that arises in the application of equation (1) is that in the problem 
treated, it is not possible to obtain an analytic expression for F(p) that 

9f 

can be differentiated to obtain — . However, if a computer program exists for 

9p 

ccanputing the solution F, then a "base" solution Fq can be calculated for a 
specified value Pq of the parameter p, and a second solution F-) = Fq + Af 
can be calculated for an incremented value p] = Pq + Ap of the parameter. 

This second solution is called a calibration solution. The derivative in equa- 

Af 

tion (1), evaluated at pQ can then be approximated by — . Consequently, 

Ap 

equation (1 ) can be used to obtain the variation of this solution from the 
base solution for any other small variation in p from its base value /if 

Af \ V 

A 2 P = P2 ~ Po then F 2 = Fq + — <^2? • Results given in reference 1 indicate 

Ap J 

that the accuracy of solutions obtained in this manner is often surprisingly 
good, even when sizable increments in the parameter are taken. 
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Methods for treating problems with discontinuities have received much 
attention in the literature. It should be noted that, since the derivative in 
equation (1 ) is taken with respect to a parameter and not with respect to the 
independent distance variable, the presence of the shock itself does not pre- 
clude a solution by equation (1 ) . It is the movement of the shock with changes 
in the parameter that causes the difficulty. Nixon (ref. 2) has treated this 
problem by using a coordinate-straining method that has been combined with 
interpolation by equation (1 ) in reference 1 . In the following development, 
discontinuous flows will not be discussed as such, under the assumption that 
for such problems the solution obtained can similarly be combined with Nixon's 
coordinate-straining method (ref. 2) . 


Theory for Solutions Depending on an Independent Function 

By extending the interpolation procedure of reference 1 with a function- 
space development, one can obtain solutions for a very general class of prob- 
lems. This class of problems involves perturbations about a base solution which 
represents a more complex flow than the simple undisturbed free-stream flow. 

Consider, for example, the nonlinear potential equation for two-dimensional 
compressible inviscid flow 


(q 2 _ 4)jj^) 4>jjy — 2(|)jj.(j)y4)jjy + (c^ 


^y^)4*yy ~ ^ 


( 2 ) 


One can assume a solution in the form of a perturbation 6(|) about a base solu- 
tion (})q 


([) = (j)o + 64> (3) 

If equation (3) is substituted into equation (2) and only linear terms in 6(j) 
are retained, then the resulting differential equation is a linear equation 
for 6<t). However, the coefficients of several terms involve either (|)q or its 
derivatives, and consequently the equation cannot generally be solved analyti- 
cally. It only becomes tractable in the case in which (IJq represents the 
constant-velocity undisturbed free-stream flow. However, perturbation solutions 
about this flow pertain only to slender configurations, and such solutions have 
already been treated extensively. 

Consider, however, problems involving nonslender configurations. Solutions 
for such problems normally are obtained by a numerical procedure. If the solu- 
tion represents the pressure or velocity distribution on a two-dimensional sur- 
face then it can be expressed as a function of the streamwise variable x. 


F = F(x) 
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The solution can be expressed as an operator Q (in general, nonlinear) operat- 
ing on the function f (x) that describes the boundary shape 

F(x) = Q[f(x)] 


Let h(x) represent a variation of the specific boundary function fo(x) . Then 
a differential of the operator can be defined with respect to function-space 
variations of f (x) . Formally, if the limit 


Q[fo + eh] - Q[fJ 

lim = L[fQ, h] 

e^O ^ 


exists in a neighborhood of if it is continuous in f at fQ, and if it 

is continuous in h at h = 0, then L is called the Gateaux differential of 
Q, and it is linear operator operating on the variation h (ref. 3, sec. 3.1). 

Although the present theory is oriented toward computer solutions, this 
concept of the differential of an operator can be illustrated by a simple ana- 
lytic example. Suppose that Q(x) were expressible in the form 

Q[f (x)l = f (X) [f (e)l 2 de 


Replacing f by fo + and taking the limit 


Q[fo + £h] - Q[fo] 
lim 


yields the result 

de + f(x),J 2f(e) h(e) de 
0 


which is linear in h. 

When Q is not expressible in analytic form, but only as the output of 
a computational procedure, then, in a manner somewhat similar to that of ref- 
erence 1 for incrementing a parameter, this operator can be estimated by 


Q = h(x) 


I 


X 


[f (e)] 
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obtaining the difference between a calibration solution and a base solution for 
slightly different boundary functions. Since L is linear, the solutions for 
other variations proportional to h can be obtained from the formula 


Q[fo + ah] = Q[fo] + L[fo,ah] = Q[fc,] + aL[fo,h] 


(4) 


This procedure, although somewhat analogous to that of reference 1 for variation 
of a parameter, differs from it in several notable respects. The first obvious 
difference is that the change in the solution is accomplished by varying a func- 
tion rather than incrementing a parameter. Secondly, the differential of Q 
is calculated directly, without, in general, computing a quantity analogous 
to the derivative in equation (1 ) . In fact, such a function-space derivative, 
called the Gateaux derivative Q’[f] = G[f], can be defined by the formula 


L[f,h] = G[f]h (5) 


However, when h is given and L is estimated by comparing base and calibra- 
tion solutions, the calculation of G[f(x)] by dividing by h is generally 
not applicable. Exceptions to this situation will be discussed in a later 
section. 

If it were possible to compute G[f], then the effect of making an arbi- 
trary variation in f could be used to obtain the effect of making any other 
type of variation. In the usual case, however, only the solution corresponding 
to boundary variations proportional to h(x) can be obtained from a single 
calibration solution. If a different type of variation is considered, then a 
second calibration calculation is required. Thus, for the two types of varia- 
tions, equation (4) gives 


Q[fQ + a-i h-| ] s ctf^] + a-|L[fQ,hi] 

Q[fo + 32h2] = Q[fol + 32L[fo,h2l 

Furthermore, since L is a linear operator on h. 


Q[fo + a-ih] + a2h2] = Qtfo^ + aiL[fQ,h-|] + a2l.[fo,h2] 


( 6 ) 


Thus, if two calibration calculations have been performed, a large variety of 
shapes determined by the various linear combinations of h-] and h 2 can be 
analyzed. Clearly this procedure can be extended to any number n of linearly 
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independent variations (i = 1 , n) provided that the total variation 

n 

^ remains small. Thus, even in the absence of analytic expressions 

i=l 

a partial synthesis of the operator L can be generated. 


Application of Theory 

Airfoil pressure distributions .- To illustrate the application of the 
method, two types of problems have been selected. The first problem involves 
the calculation, by the method of reference 4 of the upper surface pressure 
distribution on an airfoil at M = 0.56 and a = 2.0°. The base design, with 
its pressure distribution, is shown in figure 1(a). The airfoil is 17.5 per- 
cent thick with its maximum upper surface ordinate at the 30-percent-chord 
station. A variation fo + h-j of the upper surface shape is shown in fig- 
ure 1 (b) , together with its pressure distribution. This variation is 17.3 per- 
cent thick with its maximum upper surface ordinate at the 40-percent station. 

A second variation fQ + h 2 is shown in figure 1 (c) . It is 1 9 percent thick 
with its maximum upper surface ordinate at the 30-percent station. 

The first variation was obtained by shifting the abscissa, 


X = X - 0.476x(l - x) 


and defining the new surface by 
fo h] = fQ (x) 

Thus, 


h] = fo(x) - fo(x) 

The other variation was defined by 


h 2 = 0.04x(l - x) 


Mciny possible airfoil shapes can be obtained by changing the ordinate of 
the fo design in increments proportional to h-| , to h 2 » or to seme linear 
combination of the two. An example is shown in figure 2. It is represented 
by the combination fo + 0.6h-| - 0.4h2r which is 16.7 percent thick with its 
maximum upper surface ordinate at the 37-percent station. 
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Figure 2 also shows a comparison of the pressure distribution computed by 
the full nonlinear algorithm with that obtained by the present method. It is 
seen that, for this example, the accuracy of the approximation is very good. 

The actual pressure calculations for this example are given in table I. 

For many problems, the output function is not uniformly sensitive to small 
variations in the input function. For example, airfoil pressures tend to be 
sensitive to variations in the geometry very near the nose. If, for a particu- 
lar case, this sensitivity causes a problem, then smaller variations in the 
geometry shape must be taken. 

Supersonic forebody pressure distributions .- The second type of problem 
treated is that of computing the pressure distribution on pointed axisymmetric 
forebodies at supersonic speeds. Figure 3(a) shows the base configuration 
fo(x) with its pressure distribution, calculated by the method of reference 5, 
for a free-stream Mach number of 3.0. These examples are relatively straight- 
forward, and not representative of the large complex computer calculations for 
which the method is intended. They are, however, adequate for purposes of 
illustration; and they possess the further advantage that they do not involve 
specialized knowledge or notation, as would be required, for example, in a 
boundary-layer calculation. 

Figures 3(b) and 3(c) give the body shapes and pressure distributions, 
respectively, for two variations: fo + h-] (h-| = 0.04x(1 - x) ) , and £q + h 2 

(h 2 = -0.02 sin ttx) . Figure 4(a) gives results for the linear ccxnbination 
fQ + 0.3h-| + 0.7h2, both by direct calculations and by linearly combining the 
pressure distributions of figures 3(a), 3(b), and 3(c). Figure 4(b) gives a 
similar comparison for the combination fQ + 0.3h-| - 0.7h2. In both cases, 
the approximate calculations are nearly indistinguishable from the full non- 
linear results. 

The magnitude of the variation that can be taken with the calculation 
remaining in the linear range depends on the type of problem considered. If 
small increments in the input function f (x) cause large changes in the out- 
put function F(x), then the variations taken must be very small. If small 
increments in f (x) cause discontinuous jumps in F(x), then the linear theory 
is not applicable. 

An example of a sizable variation h 3 (h 3 = -0.02 sin 2t7x) of the base 
configuration fo(^) (fig. 3(a)) is shown in figure 5. The resulting pres- 
sure distribution differs radically from that of the base shape. In this 
case, the variation cannot be considered to be an incremental change in the 
base design. Nevertheless, when it is treated as such, as in the combination 
fQ + 0.3h-| + 0.7h3 shown in figure 6, the results show remarkably good agree- 
ment with the exact calculation. 


Implementation of the Method 

Application of the theory described in the previous section requires some 
simple software in addition to the basic computer algorithm that computes the 
result F(x) = Q[f(x)] from the input distribution f(x). For each case, the 
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input X and f arrays as well as the output x and F arrays must be stored 
with distinguishing case labels. Then a short auxiliary program can perform the 
following functions: 

1 . Read the x, f, and F arrays for each case. 

2. Establish a common array of the independent variable x for all cases, 
by interpolation, if necessary. 

3. Distinguish one configuration as the base fQ(x). 

4. Calculate the incremental functions L[fQ,hil by obtaining the dif- 
ferences Q[fo + h^] - Q[fol = LCfofh^] 

5. Compute and output results for any linear ccxnbination 


Q[fo ■*"^3ihi] — Q[fQ] + ^ a [ fQ, h-[] 


Considerations on Approximating Gateaux Derivative 

It was pointed out previously that the function- space differential is nor- 
mally computed directly, whereas, for an ordinary function of real variables, 
the differential is defined in terms of the derivatives (eq. (1)). The question 
arises as to whether a function-space analogue of the derivative exists. As 
mentioned earlier, under certain conditions for the existence of the linear 
operator L[fQ,h] (ref. 3, sec. 3.4), the Gateaux derivative G = Q' of the 
nonlinear operator Q can be defined by 

L[fo,h] = G[fo]h 


If L[fofh] were known in the form of an analytic expression containing a fac- 
tor of h, then an expression for G could be obtained by a simple division. 

The linear differential L could then be computed for any variation fQ + h 
by this equation. Thus, the linear approximations to the solutions for all con- 
figurations near fQ could be obtained with only two calculations. 

For the type of problem under consideration, in which the operators are 
not expressible in analytic form, such a procedure is generally not applicable. 
However, in the next section an example will be provided for which the deriva- 
tive formulation appears to be useful. 


Example of Approximating the Derivative 

The advantages of the derivative formulation of the problem (eq. (5)) pro- 
vide motivation for exploring the possibility of obtaining an estimate for the 
derivative when the solution for F is provided in the form of a numerical com- 
puter procedure. When a problem can be formulated in several ways, one of these 
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possibilities may lead more naturally to a relation of the form of ecjuation (5) 
from which the function-space derivative can be approximated. Of course, in 
such a procedure division by h at its zeros must be avoided, with values of 
the derivative at these points obtained by interpolation. 

Consider, for example, the above illustrative problem of computing the 
pressure distribution on a supersonic body of revolution. For this problem, 
the input function f (x) can be chosen to be the slope distribution of the body 
meridian line. Since it is well known that small local changes in flow angle 
are approximately proportional to the corresponding changes in pressure at super- 
sonic speeds, one might expect that the quantities h(x) and L[fQ,h] repre- 
senting the variation in body slope and the corresponding variation in pressure, 
respectively, would vary in a similar manner. Thus, a valid approximation to 
the derivative according to equation (5) appears possible. 

In order to illustrate this procedure, a calculation was performed using 
the same base body shape as in the preceding example. The variation 
-0.01 sin 2ttx was arbitrarily chosen. The slope variation is 
h-| = -0.02 cos 2itx. The result L[fo,h-|] was divided by h-] to provide 
an estimate of the derivative G[fol in accordance with equation (5). A 
second variation 0.04x(1 - x) which corresponds to a slope variation of 
h 2 = 0.04(1 - 2x) was then chosen, and the pressure increment was estimated 
by multiplying the function G by h 2 (x). Finally, for comparison, the pres- 
sure increment was determined exactly by calculating directly the pressure for 
the slope fo + h 2 . 

The results are compared in figure 7(a) for the pressure increments, and 
the total pressures are compared in figure 7(b). It is seen that the approxi- 
mate method provides a good approximation to the exact theory for this partic- 
ular problem. Such calculations with the function-space derivative should be 
useful when a rapid rough approximation is acceptable and in certain design 
problons. 


Limitations of the Theory 

It should be clear that any theory involving differentials or small vari- 
ations are limited to problems for which the output quantities are continuous 
functions of the input quantities that are to be incremented. If the output 
is relatively sensitive to variation in the input, then relatively small incre- 
ments must be taken. However, if the solution for the incremented input differs 
qualitatively fran the base solution, then the method generally fails. Scane 
examples for which the latter problen occurs are the sudden appearance of a 
shock, the appearance of leading-edge separation, and the movement of a shock 
wave or a separation point. In the case of shock movement, the linear theory 
can be applied in combination with the coordinate straining procedure of Nixon, 
which has already been noted. This procedure involves applying a shift operator 
to the incremented input function in order to "line up" the shocks for the base 
and incremented solutions. The coordinates are then shifted back a proportional 
distance to obtain the interpolated solutions. 



other Types of Problems Amenable to Treatment by Local Superposition 

Although the preceding discussion is restricted to problems for which a 
change in the output function is effected by varying an input function, a simi- 
lar treatment is applicable to other kinds of problems. In reference 1 , a 
change in the output function resulted from incrementing an input parameter, 
such as Mach number, angle of attack, or thickness ratio. Using the present 
method, it would be possible to compute the effect of incrementing each of these 
parameters separately, and then approximate the effects of incrementing various 
combinations of them simultaneously by linearly superimposing the separate 
increments. 

Another type of problem involves the increment in a single output quantity 
as a result of incrementing an input parameter. Examples of such problems 
include computing the lift of an airfoil or the wave drag of a supersonic con- 
figuration as a function of Mach number, angle of attack, or thickness ratio. 
Since, in this case, the output quantity is a simple function of several vari- 
ables, the partial derivatives with respect to each independent variable can be 
approximated, and the total differential due to varying combinations of them 
can be ccmputed by the linear nature of the derivatives. 

Suppose, on the other hand, that the increments in an output quantity such 
as lift or drag were due, not to incrementing a parameter, but to varying the 
configuration shape. The output quantity is then denoted a functional (ref. 3, 
p. 9), depending on the input function that describes the boundary shape. By 
computing the increments due to variation in the boundary shape, one can obtain 
the effects of linearly combining the variations. An interesting problem of 
this type is that of studying the changes in the effective perceived noise level 
at a ground point due to variations in the flight pattern or engine operation 
of a jet airplane. 

Finally, it should be pointed out that the inputs to this type of procedure 
are by no means limited to those obtainable from an analytic formula, or a compu- 
tational procedure, but could just as well represent experimental data. 


CONCLUDING REMARKS 

A theory has been presented for reducing the usage of large computational 
algorithms. Such algorithms are treated in this theory as nonlinear operators 
operating on an input function. The method requires a procedure for obtaining 
computer solutions for a base problem, and for a limited number of variations 
of this problem. The solutions of other problems that differ slightly from the 
base problan are then found by forming linear combinations of the increments 
obtained for the computed variations. The examples given involved nonlinear 
analysis of subsonic airfoils and of supersonic forebody shapes. 

For a restricted class of problems, a function-space derivative can be 
approximated. For such cases, only one variation from the base problem is 
required in order to approximate solutions for general variations. The example 
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that was described treated the solution for the pressure distribution on an 
axisymmetric supersonic forebody as an operator on the slope distribution of 
the body meridian. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 31 , 7 981 
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TABLE I.- PRESSURE CALCULATIONS FOR AIRFOIL DEFINED BY Fq + 0.6hi - 0.4h2 


Station, x 

Cp,o 

Cp,l 

Cp,2 

^Cp 



^2Cp 

Cp + 0.6 A-]Cp - 0.4 A2Cp 

Exact 

0.05 

-1 .130 

-0.955 

-1 .140 

0.1 75 

-0.01 0 

-1 .021 

1 

• 

o 



.10 

-1 .400 

-1 . 01 5 

-1 .566 

.385 

-.1 66 

-1 .103 

-1 . 090 

.14 

-1 .406 

-1.033 

-1 .674 

.373 

-.268 

-1 .075 

-1 .088 


-1.285 

-1 .051 

-1 .565 

.234 

-.280 

-1 .033 

-1 .055 


-1 .1 64 

-1 .0575 

-1 . 41 0 

.1065 

-.246 

-1 .002 

-1 . 01 6 


-1 .038 

-1 . 051 

-1 .250 

-.013 

-.21 2 

-.961 

-.968 


-.916 

o 

U) 

-1 .1 05 

-.115 

-.189 

-.909 

-.913 


-.795 

-.997 

-.970 

-.202 

-.175 

-.841 

-.847 

.50 

! m 

-.880 

-.730 

-.273 

-.1 23 

-.722 

-.717 

.60 


-.735 

-.530 

-.295 

-.090 

-.581 

-.570 

.69 


-.570 

-.365 

-.269 

-.064 

-.437 

-.430 

.80 

-.135 

-.339 

-.1 67 

-.204 

-.032 

-.245 

-.240 

.90 


-.069 

.043 

-.107 

.005 

-.024 

-.020 
























(b) First variation h-| (x) . 
Figure 3.- Continued. 
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Exact 




(a) First composite body £q + 0.3h-] + 0.7h2. 

Figure 4.- Comparison of exact and approximate theories for two com- 
posite configurations synthesized from the base configuration 
fg and the two variations h] and h 2 « M = 3.0; a = 0°. 
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Figure 6.- Comparison of exact and approximate theories for the 
composite configuration fQ + O.Sh-j + 0.7h3. M = 3.0; 
a = 0°. 
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